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Abstract 



We consider continuous observation of the nonlinear dynamics of single 
atom trapped in an optical cavity by a standing wave with intensity modu- 
lation. The motion of the atom changes the phase of the field which is then 
monitored by homodyne detection of the output field. We show that the 
conditional Hilbert space dynamics of this system, subject to measurement 
induced perturbations, depends strongly on whether the corresponding clas- 
sical dynamics is regular or chaotic. If the classical dynamics is chaotic the 
distribution of conditional Hilbert space vectors corresponding to different 
observation records tends to be orthogonal. This is a characteristic feature of 
hypersensitivity to perturbation for quantum chaotic systems. 
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I. INTRODUCTION 



The study of quantum nonlinear dynamics, especially systems which classically exhibit 
Hamiltonian chaos, has recently begun to focus on the response of such systems to external 
sources of noise and decoherence This direction was prompted by the observation that 
nonintegrable classical systems, when quantized will exhibit dynamics that departs from 
that expected classically on a very short time scale so short that even macroscopic 

systems should show observable quantum features in their motion . There have now been 
numerous experimental observations of the short time deviations between quantum and 
classical dynamics of nonlinear systems ||7|j8|j9|,p!0|l . The nonlinear dynamics of cold atoms in 
optical dipole potentials has proved to be a particularly fertile field for quantum nonlinear 
dynamics. Recently the effects of decoherence in quantum chaotic dynamics was studied 



using cold atoms [|Tl|,|T3)|Il] • However in all experimental observations so far, the results 
were obtained from an ensemble of systems, not from repeated observations on a single 



quantum system. Recent progress in single atom dynamics in small optical cavities 14 



indicate that it will soon be possible to study the quantum nonlinear dynamics of a single 



quantum system subject to repeated measurements [|T5| and it is towards describing such 
systems that this paper is directed. 

It is in the context of such single system dynamics that the information approach of 



Schack and Caves [0 based on hypersensitivity to perturbation becomes significant. In that 
approach the response of classical and quantum nonlinear systems to external perturbations 
is considered. In particular they show that for a chaotic system it requires a huge amount 
algorithmic information to track the classical (phase space trajectory) or quantum (Hilbert 
space vector) of a single chaotic system when it is subjected to small external perturbations. 
It is better in such cases to average over the perturbation and pay a much smaller cost 
in von Neumann entropy. In the quantum case the signature of this hypersensitivity to 
perturbation has been shown to be the distribution of Hilbert space vectors resulting from 
dynamical sequences with different perturbation histories. While this approach seems to offer 
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considerable insight into quantum and classical chaos, it is far from clear what it means for 
an experiment where the dominant source of perturbation is likely to be the measurement 
back action associated with the attempt to continuously monitor the dynamics. 



In reference an attempt was made to study hypersensitivity to perturbation arising 
from quantum measurements made on a single quantum system: a nonlinear kicked top. In 
that study the measurements were not continuous in time but rather a sequence of discrete 
readouts applied at the same time as the kicks. The results confirmed in general terms 
the observation of Schack and Caves for measurement induced perturbation. Specifically 
it was shown that if a system was initially localized on a chaotic region of phase space, 
the Hilbert space vectors resulting from different measurement histories tended to become 
orthogonal, while for initial regular states the Hilbert space vectors for different histories 
tended to remain closer together. In this paper we extend that study to the case of a 
continuously monitored single quantum nonlinear system: a single atom trapped by an 
intracavity optical dipole field. The motion of the atom changes the phase of the cavity field 
which may be monitored using phase sensitive detection of the light leaving the cavity. Chaos 
is introduced by externally modulating the intensity of the light inside the cavity. We use 



the now established techniques of quantum trajectories |jIB| , P^ to study the distribution of 



Hilbert space vectors for different measurement histories. Previous studies that use quantum 
trajectories to describe the dynamics of open quantum nonlinear systems include the work 



of Brun et al. 20 



We use the nonlinear stochastic Schrodinger equation to parallel the discussion in ref- 



erence 



p!7[] based on how information is extracted. The nonlinear Schrodinger equation 
describes a true measurement in which actual information about the quantum state of the 
monitored system is extracted from the external field. Our results confirm that a chaotic 
system, subject to different continuous observation histories, will produce a distribution of 
states that tend to be orthogonal. This means that a very tiny error in recording the mea- 
surement history will suggest a final state that is very likely orthogonal to the actual final 
state. In this way the intuitive idea that chaos constrains predictability is carried over to 



continuously observed single nonlinear quantum systems. 



II. THEORETICAL MODEL OF HOMODYNE MEASUREMENTS ON SINGLE 

ATOM DYNAMICS IN CAVITY QED 

Recently, experiments in cavity QED have achieved the exceptional circumstance of 
strong coupling, for which single quanta can impact the atom-cavity system. The trapping 



of a single atom in a high- finesse cavity has been realized [jl4[. In these systems we must 
treat quantum mechanically both the optical and electronic degrees of freedom as well as the 
center-of-mass motion of the atom. In our model the atom is in the optical dipole potential 
of a cavity standing wave which is blue detuned from an atomic resonance so that there is 
a net conservative force acting on the atom in the direction of decreasing intensity. This 
interaction does not change the intensity of the optical field but it does change the phase 
by an amount that depends on the atomic position. As the atom moves in the cavity it 
changes the phase of the field and if this phase change can be monitored we can effectively 
monitor the atomic position. This can be accomplished by a homodyne measurement of 
the field leaving the optical cavity. Mabuchi et al ||23| have already demonstrated this kind 



of measurement at the level of a single atom. A similar model for an atom trapped in a 



harmonic optical potential was recently discussed by Doherty et al 

The basic theoretical description can be given as master equation for a two-level atom 
coupling to a single electromagnetic mode via the Jaynes-Cummings interaction Hamilto- 



nian, including the quantization of the atomic center-of-mass pi| , |22| . The Hamiltonian in 
Schrodinger picture can be written as 

H = ^ + hujA<y^(T' + huca'^a + hEo{ae~''^''* + a+e*"^-^*) + kg sm{kLx){aa^ + a+cr"), (1) 

where p is the momentum of atom, M is its mass, uja, ^c-, and ujl are the the two-level 
resonance frequency, the cavity frequency, and the frequency of the driving laser field, re- 
spectively. The term Eq is a constant proportional to the amplitude of the driving field, g 
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is the coupling constant of the interaction between driving field and atom, cj^ and a" are 
the raising and lowering operators for the two-level atom, and and a are the creation 
and annihilation operators for the cavity field. We assume that the detuning A is positive 
and A = — Ci^L ^ fif, r and Uc = ujl, where T is the atomic dipole decay rate. In the 
interaction picture the Hamiltonian can be simplified as 

H' = Heff + hEo{a + a+) (2) 

where 

Heff = — a"'"a cos(2A;ra;), (3) 



is the effective Hamiltonian ||21|. Note that the effective interaction does not include the 
driving laser field. 

We denote A the density operator for the joint state of the atom and the cavity. Then 



the master equation for A is p^,|25 



dK 1 - K 

— = —[Heff, A] - iEo[a + a+. A] + -(2aAa+ - a+aA - Aa+a), (4) 

where k is the cavity decay rate. Note that if the cavity is driven by a strong coherent 
field and if it is strongly damped at the rate k, the field state will relax to approximately a 
coherent state with amplitude a = — 

We assume that Eq/ n <^ 1. Then we can transform the total state by 

A = D+{a)KD{a). (5) 



Therefore a ^ a + a and a+ ^ a+ + a*. We then expand A |24 



A = po ® |0)a(0| + (pi ® |l)a(0| + H.C.) +P2® |l)a(l| + {p 2 ® |2)a(0| + H.C.). (6) 

The reduced density operator is p = Tr(A) = po + P2 and the master equation after 
adiabatic elimination is 

|._Lp,,|_fl|J,li,p]l. (7) 



Here 



= ^ (8) 



is the diffusion constant and 



J = - cos{2kLx), (9) 



"2 

Ho = ^ + nxl (10) 



where 



X^^. (11) 
The conditional master equation for the optical field undergoing continuous Homodyne 



measurement is p^ , |25[] 



,dpc X I I I X /—dWit) , \ /ir.\ 

(— ) = -[2apca^ - a^apc- PcO^a) + ^/k — -— (apc + Pc^ - (« + « )cPc), (12) 
at field I at 

where dW(t) is the infinitesimal Wiener increment. In this equation pc is the density 
matrix that is conditioned on a particular realization of the Homodyne current up to time 
t. The corresponding stochastic Schrodinger equation is 

d\^,{t)) = dt[~iH,ff/h - ]^Ka+a + I{t)a]\^,{t)) , (13) 

where 

/(t) = «:(a + a+) + v^^^ (14) 

is the measured current. Using Eq. (7), we can derive the nonlinear stochastic Schrodinger 
equation by adiabatic elimination, 

d\i,,{t)) = dt[-tHo/h - DJ^ + lA{t)J]mt)), (15) 
where /a = 4D(J) + V2D^^ 



dt 



The normalized nonlinear stochastic Schrodinger equation is 

d\Mt)) = dt[-zHo/h - D{J - (J),)2 + v/2D(J - (J),)^^]|^,(t)). (16) 

Given the modulation frequency uj, we can define dimensionless parameters by t = ut, 
P = il^)P^ ^ = 2fcLa;, ^0 = wi^Ho, g = g/uj, E = E/uj, A = A/a;, R = k/uj, D = D/uj, 
and X = x/^- This yields the commutator relation 

[x,p]=ik, (17) 

where k= is the dimensionless Plank constant. 

Omitting all the tildes , the equivalent equations are similar except h is replaced by k 
and the dimensionless Hamiltonian is 

Ho = — -^cosx, (18) 

where ^ = ^giT^X- Iii order to study chaos in a quantum system, we consider a periodic 
modulation of the driving field Eoit) = Eo\^l — 2e cost. The expressions of the stochastic 
equations (14) and (15) will not change except that h is replaced by k and D replaced by 
D{1 — 2ecost) and ^ in Eq. (17) is replaced by ^(1 — 2ecost). 



III. SENSITIVITY TO DIFFUSION CONSTANT 

We assume that initially the atomic center-of-mass wave function is in a Gaussian mini- 
mum uncertainty state with the position representation 

We take Xq = 0, po = 1.0 as for these values the state is localized on a second order 
period one resonance and is thus localized in a regular region of phase space (see Fig. 1). 



For ax = 0.3906, k = 0.25, ^ = 1.2, Dyrting et al. ||26| have shown that the system will 
coherently tunnel between the two corresponding second order period one resonances. We 
use a Split Operator Method and FFT (Fast Fourier Transformation) |2^ to obtain 



the numerical solution of the stochastic Schrodinger equations. In this scheme the kinetic 
operator and potential operator are used separately to propagate the wave function: 

exp[-imt/k] ~ exp[-t{Py6t/4k] exp[-i(V)St/2k] exp[-t{Py6t/4k]. (20) 

The computing errors are of 0{6t^). Here V is the effective potential which includes a 
stochastic term, 

^2n , dW{t)\ /;77^,j ,dWit)^ 



V = -at) cosi - tk[D{t){J - (J),)^(l + - pD{t){J - {J)c)^^]. (21) 

where ^(t) = ^(1 - 2ecost) and D{t) = D{1 - 2ecost). We introduced the ^^^^^ term to 
keep the expression consistent with the normalized nonlinear stochastic Schrodinger equation 
after an expansion of the exponential function. 

In order to compare the quantum and semi-classical stochastic evolutions, we calculate 



the Wigner function [|^, |^ 

1 /■ J / y, , , y 



If y y 

^(a;,p) = ^y dy{x - -\p\x + -) exp{ipy/k). (22) 



This expression can be interpreted as the Weyl- Wigner correspondence of the density 
operator. To give the dynamical equation for the Wigner function that is the quantum 
correspondence of a classical Liouville equation we use the Weyl- Wigner correspondence of 



an operator F = AB which is pO,31 



F(x,p) = A{x,p)XB{x,p), (23) 

where X = exp[^(^^ — ^^)] and the arrows on the operators denote the term on which 
the operator is to be applied. Alternatively, we obtain 

F(x,p) = A(x-||,p+||)S(x,p). (24) 

When we apply this formula to the products appearing in the Master equation we can 
readily obtain the phase space equation for the Wigner function we are looking for 

dP ^ djmdp _ djmdp ^ j,^.^,^2^d^^ .25) 

dt dq dp dp dq dp"^ ' 
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where -ffo(^) is the classical Hamiltonian including modulation HQ{t) = ^ — ^(t) cosx. We 
give the corresponding classical stochastic F-P equation from quantum nonlinear stochastic 
Schrodinger equation |^ 

^ = -^(t) sinx + ^2D{t)ksmx ^^J'^\ (27) 



To describe the classical distribution we use the classical Q function The initial state is 



a bivariate Gaussian centered on (xo,Po) with position variance 6x and momentum variance 



Qo ix,p) = - — j== exp — exp — , 28 



h2 h2 

where the classical variances 6^ and 6p are related with quantum parameters = ^ + 

6p = — h Cp. The evolution of Q function is Q{x,p,t) = Qo[x{x,p, —t),p{x,p, —t)], where 

x{x,p, —t),p{x,p, — t) is the trajectory generated by Hamilton's equations. 

To compare the quantum dynamics with the classical conditional dynamics, for quantum 
system we study the ensemble with the same initial condition but with random trajectories. 
It shows that when D is very small, for Homodyne measurement, the evolution of average 
momentum (p) and average variance of momentum ((p^) — (p)^) for the ensemble show 
coherent tunneling. Therefore the perturbation is not serious for small D when the initial 
state is in the regular region of the classical phase space. 



For the classical dynamics for small D the results are close to the no diffusion case ^6 
Obviously we therefore expect that we obtain different results between classical and quantum 
conditional dynamics (Fig. 2). 

However, if the diffusion constant is large enough, we obtain almost the same result as in 
the classical case (Fig. 3). For a single stochastic measurement, the terms in the normalized 
nonlinear Schrodinger stochastic equation due to the measurement depend on the quantity 
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(J — {J)c)- We expect that for some range of values of D, the stochastic measurement terms 
would drive the system towards an oscillating trajectory for which 

{J\ ^ {J)l (29) 

Therefore for the ensemble which includes many random trajectories, the results will ap- 
proach that of the classical conditional dynamics. 



IV. SENSITIVITY TO CHAOTIC AND REGULAR INITIAL STATES 

We again assume that the wave function is initially in a minimum uncertainty state in 
the position representation. We choose two initial locations in classical phase space. In the 
first case, xq = 0, = 1-0 it is in the regular region of classical phase space. In the second 
case xo = —2.5, po = 1.0, it is in the chaotic region (see Fig. (1)). 

Because the measurement will perturb the quantum state, we hope to be able to compare 
the effect of measurement noise of different trajectories. Here the angle 6ij{t) is defined 
between two normalized state vectors I'lpiit)) and li'jit)) as |[16|jr7[ 



9ij{t) = cos-^ \{ij,{t)\^,{t))\. (30) 
In the position representation, this is 

liUtMjm = I / ^iix,tyi;,{x,t)dx\. (31) 



As a measure of the distribution of the state vectors in Hilbert space we can calculate the 
average angle between all pairs of vectors. We define 

^aUt)= ^^,^_^^ T.O^At)^ (32) 

where N is the number of trajectories. 

In Fig. 4, we plot Oaveif) for the two initial states mentioned above, evolved up to 200 
cycles. We used up to 40000 steps and N = 1000 trajectories for our calculation. As can be 
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seen, the average angle between vectors starting in the chaotic phase space region is larger 
than that of the regular initial state. 

If the system started in the regular region, small errors in recording the results of the 
measurement will not be a very serious problem because the conditional states form trajec- 
tories which will remain close in Hilbert space. Therefore, if we consider the distribution 
of Hilbert angles at a fixed strobe number (Fig. 5), the distribution is centered at a small 
angle for an initial regular state. On the other hand, for an initial chaotic state, the peak 
location approaches 71/2 (Fig. 6). It means that for an initial chaotic state most vectors are 
far apart from each other. Therefore if the initial state was in the chaotic region of phase 
space, it will be much more difficult to infer the system state reliably from the measurement 



results. This result is consistent with the results of the quantum kicked top [|T7 . 

In summary, we have demonstrated that for continuous Homodyne measurement of sig- 
nals from the quantum system of single atom dynamics in cavity QED, the measured results 
are influenced by the diffusion constant and whether the initial states are in regular or 
chaotic phase space regions. We have shown that if the diffusion constant is large enough, 
the average measurement results are similar to classical conditional dynamics and for small 
diffusion constant the initial chaotic state will be more sensitive to errors in recording the 
measurement results than the initial regular state. 

V. ACKNOWLEDGMENT 

One of authors(XML) would like to thank Dr. J. F. Corney and Dr. M. Gagen for useful 
discussions on conditional master equations and Homodyne measurements. 



11 



REFERENCES 

[1] W.H.Zurek and J.P.Paz, Phys. Rev. Lett. 72, 2508 (1994). 

[2] R.Schack, G.M.D'ariano and C.M.Caves, Phys. Rev E 50, 972 (1994). 

[3] S.Fishman, D.R.Gremple and R.E. Prange, Phys. Rev. Lett. 49, 509 (1982). 

[4] D.L.Shepelyansky Phys. Rev. Lett. 56, 677 (1986). 

[5] B.V.Chirikov, F.M.Israilev, D.L.Shepelyansky, Physica D 33, 77 (1989). 

[6] W.H.Zurek, Physica Scripta T76, 186, (1998). 

[7] E.J.Galvez, B.E.Sauer, L.Moorman, P.M.Koch and D. Richards, Phys. Rev. Lett. 61, 
2011 (1988). 

[8] J.E.Bayfield, G.Casati, I.Guarneri and D.W.Sokol, Phys. Rev. Lett. 63, 364 (1989). 

[9] M.Arndt, A. Buchleitner, R.N. Mantegna and H.Walther, Phys. Rev. Lett. 67, 2435 
(1991). 

[10] F.L. Moore, J.C. Robinson, C.F. Bharucha, Bala Sundaram, and M.G. Raizen, Physical 
Rev. Lett. 75, 4598, (1995). 

[11] B. G. Klappauf, W. H. Oskay, D. A. Steck, and M. G. Raizen, Physical Rev. Lett. 81, 
1203, (1995). 

[12] B. G. Klappauf, W. H. Oskay, D. A. Steck, and M. G. Raizen, Physical Rev. Lett. 81, 
4044, (1995). 

[13] H. Ammann, R.Gray, I Shvarchuck, and N. Christensen. Physical Rev. Lett. 80, 4111 
(1998). 

[14] J.Ye, D.W.Vernooy, and H.J.Kimble, Phys. Rev. Lett. 83, 4987(1999). 

[15] A.C.Doherty, S.M.Tan, A.S.Parkins and D.F.Walls, Phys. Rev A 60, 2380 (1999). 

12 



[16] Riidiger Schack, and Carlton M. Caves, Phys. Rev. A 53, 3257(1996). 
[17] J.K. Berslin, and G. J. Milburn, Phys. Rev. A 59, 1781(1999). 

[18] H. J.Carmichael, An Open Systems Approach to Quantum Optics, Lecture Notes in 
Physics: New Series m. Monographs, ml8, (Springer, Berhn, 1993). 

[19] K. Molmer, Y.Castin, Quantum Semiclass. Opt. 8, 49 (1996). 

[20] T.A.Brun, Gisin N, P.F.O'Mahony , M. Rigo, Phys. Lett A 229, 267(1997). 

[21] A.M.Herkommer, H. J.Carmichael and W.P.Schleich, Quantum Semiclass. Opt. 8, 
189(1996). 

[22] Ralf Quadt, Matthew CoUett, and Dan F. Walls, Phys. Rev. Lett. 74, 351(1999). 
[23] H. Mabuchi, J.Ye, H.J.Kimble, Apphed Physics B-Lasers and Optics 68,1095,(1999). 
[24] H. M. Wiseman, and G. J. Milburn, Phys. Rev. A 47, 642(1993). 
[25] J. F. Corney, and G. J. Milburn, Phys. Rev. A 58, 2399(1998). 
[26] S. Dyrting, and G. J. Milburn, Phys. Rev. A 51, 3136(1995). 
[27] Ronnie Kosloff, J. Phys. Chem. 92, 2087(1988). 

[28] WiUiam H Press et al (1992), Numerical Recipes in C , Cambridge, University Press. 

[29] E. Wigner, Phys. Rev. 40, 749(1932). 

[30] H. J. Groenewold, Physica 12, 405(1946). 

[31] L. Wang and R. F. O'Connell, Found. Phys. 18, 1023(1988). 

[32] D. F. Walls, G. J. Milburn(1994), Quantum Optics, Springer- Verlag Berlin Heidelberg. 
[33] Wenyu Chen, S. Dyrting, and G. J. Milburn, Aust. J. Phys. 49, 777(1996). 



13 



FIGURES 

2 

FIG. 1. Stroboscopic portrait of the system with Hamiltonian Hq = ^ — C(l — 2e cost) cos x, 
where e = 0.2, ^ = 1.2. 

FIG. 2. Evolution of average momentum {p) and and average variance of momentum (jp) — (p)^ 
for classical and quantum conditional dynamics when D = 0.001. 1000 random trajectories are 
taken. Solid line, classical conditional dynamics, dashed line, quantum conditional dynamics. 

FIG. 3. Evolution of average momentum {p) and and average variance of momentum (p^) — (p)^ 
for classical and quantum conditional dynamics and D = 0.1. 1000 random trajectories are taken. 
Solid line, classical conditional dynamics, dashed line, quantum conditional dynamics. 

FIG. 4. Evolution of average angles in Hilbert space for D = 0.001. Solid line, in chaotic region 
initially xq = —2.5, po = 1.0. Dashed line, in regular region initially xq = 0.0, po = 1.0. 

FIG. 5. Distribution of angles in Hilbert space at strobe number 200 for initial regular state and 
D = 0.001. 

FIG. 6. Distribution of angles in Hilbert space at strobe number 200 for initial chaotic state 
and D = 0.001. 
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